Incommensurate superfluidity of bosons in a double-well optical lattice 



OO 
O 

o 

(N 

CD 
00 

(N 



5-H 

CD 



I 

C 

O 

o 



> 

m 

O 

OO 

o 



X 



Vladimir M. Stojanovic, 1 Congjun Wu, 2 W. Vincent Liu, 3 and S. Das Sarma 4 

1 Department of Physics, Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, USA 
Department of Physics, University of California, San Diego, California 92093, USA 
^Department of Physics and Astronomy, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, USA 
* Condensed Matter Theory Center, Department of Physics, 
University of Maryland, College Park, Maryland 20742, USA 
(Dated: September 22, 2008) 

We study bosons in the first excited Bloch band of a double-well optical lattice, recently realized 
at N1ST. By calculating the relevant parameters from a realistic nonseparable lattice potential, we 
find that in the most favorable cases the boson lifetime in the first excited band can be several 
orders of magnitude longer than the typical nearest-neighbor tunnelling timescales, in contrast to 
that of a simple single-well lattice. In addition, for sufficiently small lattice depths the excited band 
has minima at nonzero momenta incommensurate with the lattice period, which opens a possibility 
to realize an exotic superfluid state that spontaneously breaks the time-reversal, rotational, and 
translational symmetries. We discuss possible experimental signatures of this novel state. 



Optical lattices provide an exquisite tool for con- 
trolled exploration of novel types of order in cold atomic 
gases [1]. In particular, realization of a (quasi) two- 
dimensional (2D) double-well (DW) optical lattice at 
NIST was one of the latest major developments in the 
experimental cold-atom physics 2], motivating further 
experimental [3| and theoretical [4j efforts. 

Ultracold atoms, either fermionic or bosonic, in higher 
Bloch bands have recentl y ig nited a great deal of inter- 
est d H 0, i, S, M, EH, E El- In the fermion case 
Pauli blocking enables one to populate high Bloch bands 
by simply increasing the atomic density [8|. By contrast, 
the true ground-state condensation of bosons only occurs 
in the lowest s-orbital band even for high boson densi- 
ties. When the majority of bosons populates a higher 
band, bosons are in excited states with a finite lifetime. 
Isacsson et al. o| have studied bosons in the first excited 
band of ID optical lattice and found lifetimes that are 
on the order of 10 — 100 times longer than the typical 
nearest-neighbor tunnelling time, a prediction that has 
recentl y b een experimentally corroborated by the Mainz 



group 11]. In this Letter, we study weakly-interacting 



cold bosons populating the first excited band of a quasi- 
2D DW optical lattice. 

We show that in the superfluid regime Bose-Einstein 
condensation (BEC) takes place at an incommensurate 
nonzero momentum, which spontaneously breaks time- 
reversal, rotational, and translational symmetries. We 
further demonstrate that, due to vastly reduced available 
phase space for the decay to the lowest band, the lifetime 
of a Bose gas in the first excited band of a DW lattice 
can be several orders longer than the inverse tunnelling 
rate, which in turn sets the characteristic time needed to 
establish phase coherence in the system jl4j |. 

The DW lattice consists of a nonseparable optical po- 
tential in the x-y plane and a conventional optical poten- 
tial in the z-direction. In the case with DWs oriented in 
the x-direction, in a coordinate system with the origin at 



a maximum point of the "in-plane" -lattice light intensity, 
the optical potential in the x-y plane is given by Q 



V(x, y) = 2V {[cos(2fc L2 /)-cos(2fc L x)] 
+ 2r[cos(fc L x) + cos(fc L y)] 2 } . 



(1) 



Here fc L = 2ir/\ is the magnitude of the laser wave- vector, 
Vo = — | Vb| < (red-detuned lattice), while r = I z /I xy 
stands for the relative intensity of two components of 
light with the in-plane and out-of-plane polarizations. 
The band structure of the single-particle Hamiltonian 
H Q = -(ti 2 /2m b )(dl + dy) + V(x,y) (mf,-boson mass) 
in the x-y plane can be solved by using the plane-wave 
basis. The corresponding matrix elements of Hq read 
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k + G m ,„) = E R {[(k x /k h )+m + : 
= (k\H \k+G 0l±1 ) = 2rV , 



G±i,±i) = (r - l)Vb , 
G±i, T i) = (r + l)Vo, 



(2) 



where k = (k x , k y ) is the wave- vector, G m; „ = mhi+nh2 
(m, n-integers) are the reciprocal-lattice vectors (with 
basis bi, 2 = k L (e x ± e y )), and E R — fi 2 ^/(2md is the 
recoil energy. 

The dispersion of the first excited band is depicted 
in Fig. [H where the energies are expressed in units of 
E B . While the lowest band has minimum at k = 0, the 
first excited band has maximum at k = and minima 
for k 7^ 0. For larger values of |Vb|, these minima oc- 
cur at commensurate wave- vectors of k = (±l,0)fc L and 
(0, ±l)fc L . However, for optical lattice depths smaller (in 
absolute value) than some r-dependent threshold value, 
i.e., for | Vol < Vthtj), these minima occur at wave- 
vectors K that are incommensurate with the lattice pe- 
riod in both x- and y-directions and independent of V"o : 
K x = ±0.5 fc L , K y — ±0.5 For example, the threshold 
value for r = 0.08 is V th « 1.08 E R , while for r = 0.15 
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FIG. 1: Excited-band dispersion for V /E R = -1.0, r = 0.08. 




( in units of k L ) 



FIG. 2: Contour plot of the excited-band dispersion for 
Vo/Er — —1.0, r — 0.08. Band-minima are marked by 'x'. 



it is Vth ~ 0.83 E R . Contour plot of the excited-band 
dispersion for Vo/E R = —1.0 and r = 0.08 is displayed in 
Fig. [2] The band-minima are only two-fold degenerate, 
related by a mirror symmetry, because (±^fc L , ±|fc L ) are 
equivalent to each other and so are (i^&L? ~F 2^l) - 

The two relevant bands originate from the lowest- 
energy states of a DW: the lowest (+) band from the 
DW ground state (even parity), the excited (-) band 
from the first excited state (odd parity) of a DW. These 
two states can be sought in the form of "bonding" and 
"anti-bonding" linear combinations of single-well Gaus- 
sians, respectively, similar to the Heitler-London varia- 
tional approach to the molecular ion. For the DW 
comprising potential- minima at (a;i,0) and (#2,0), 



$±(x,y) 



(p(x - x x ,y) ± ip(x - x 2 ,y) 



(3) 



where (p(x,y) = (ttct 2 ) 1/2 e ^ 2 is a 2D Gaussian 
and S = J (p*(x - xx,y)<p(x - x 2 ,y) d 2 r = e~ (r /( 4<T ' 
is the overlap integral of two such Gaussians, a distance 
b = xi — x\ apart. Optimal value a of the Gaussian- 
width a is obtained by minimizing the expectation value 



($+(x, y)\H \$ + (x, y)) of the DW ground-state energy 
over this parameter. The value of <r becomes larger with 
decreasing lattice depth (i.e., for decreasing \ Vq\). In par- 
ticular, our calculation shows that for |Vo|/-E R < 1.0 one 
has <7 > 0.39 b. Strictly speaking, the Wannier func- 
tions are well described by the ansatz in Eq. only 
for not-too-small |Vo| (implying that the mid-barrier be- 
tween single wells is sufficiently high). By comparing the 
variational ground state energies of a DW corresponding 
to different lattice depths with band-structure calcula- 
tions, we can estimate that this approach is accurate for 
\Vq\ > 0AE R , which puts the lower bound on the lattice 
depths that will hereafter be discussed. 

The z-dependent part of the full 3D Wannier functions 
$±(x, y,z) — $±(x,y)(t>(z) takes on the standard Gaus- 

sian form <f>(z) = (7i"£ 2 ) _1/4 e : where £ z is the effective 
harmonic "zero-point" length in the z-direction, charac- 
terizing the transverse confinement of the system. 

The intra-band and inter-band Hubbard interaction 
parameters are given by U± = (g/^zV^n) J §±(x,y)d 2 r 
and U + - = (g/£zV2n) J $+(x, y)$_ (x,y)d 2 r, respec- 
tively, where g = 47rS 2 a s /m(,. Calculation yields 



U+- = 2(a a /^) 

E R ^{o a KY ' 

U± _ 1 + 3e- fc2 / 2CT o ± 4e- 3b2 / 8CT o 

~ (l±e^ b2 / 4CT o) 2 



(4) 
(5) 



These interaction energies are proportional to the ratio 
a s /£ z , the realistic values of which can be estimated to be 
between 0.03 and 0.07 0. Evaluation of Hubbard U% 
displayed in Fig. [31 shows that they are smaller than 
the bandgap between the lowest and the excited band 
f° r l^o I i5 2.2 E R , the latter being the largest optical 
lattice depth that we shall be concerned with in what 
follows. [It is useful to note, however, that parameter Vo 
here does not have entirely the same meaning as lattice 
depth in the case of conventional optical lattices.] For 
instance, for V /E R = -1.0, r = 0.08, and a s /^ z = 0.06, 
values of Hubbard parameters arc U+ = 0.0426 E^, U- = 
0.0498 E R , U+- = 0.0420 E R . Unlike the situation in 
ordinary (single- well) optical lattices, where the Hubbard 
energy of the excited (p) band is smaller than that of the 
lowest (s) band, here we find that U- > U+. 

While the '—'-band minima are degenerate, condensate 
fragmentation 15| in momentum space is prevented by 
the presence of interactions, no matter how weak. It is 
also easy to check that condensation into any state that 
is a linear superposition of two different band minima 
has higher energy than in a sing le minimum. Therefore, 
our state is a "simple" BEC [16| at wave-vector denoted 
with K. It is interesting to note that '—'-band has nearly 
flat dispersion around the minima, which may pose a 
challenge to the experimental realization of this state. 

Our primary objective is to study the influence of 
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FIG. 3: Hubbard interaction parameters (in units of E R ) for 
r — 0.08, a s /£ z — 0.06, and different lattice depths. 

interactions on the lifetime of bosons in a band mini- 
mum of the excited band. Thus it is essential to eluci- 
date the dominant decay process involved. It is worth- 
while noting that the DW-type systems have asymmetric 
level-spacings, with energy levels appearing in pairs and 
the energy gap between the '—'-band and the second- 
excited band (W32) being much greater than the one be- 
tween the '+'- and '—'-bands (W21): f° r instance, for 
Vq/E r = —1.0 and r = 0.08 our band-structure calcula- 
tion yields W 32 = 1.52 E R , W 2 i = 0.079 J5 B . Therefore, 
promotion of particles from the '— ' band to the yet higher 
one, relevant for the decay channel studied in Ref. |5j, is 
energetically very costly. Thus the dominant decay pro- 
cess we need consider is the one where two bosons decay 
from the '—'-band- minimum to the '+'-band. The in- 
crease in the interaction energy, which is predominantly 
due to the exchange-energy contribution (inter-band in- 
teractions, here described by U-\ ), compensates for the 

decrease in the single-particle band energy. 

We recall the expression for the interaction energy 

,(6) 

m;k,k',q 



2iVo 



where m, n, ml ', n' are band indices, and No is the number 
of unit cells in the system. For the first two bands, inter- 
action parameters are U++ = [/+, U = U- and 

U++,— = U— t++ = [/+-',_+ = u+-,+J= £/_+,_+ = 

U ^ = U-\ Let us assume that the initial state is 

= [a_ K ) Ar |vac), where N is the number of particles 
in the condensate and |vac) stands for the boson vacuum, 
while the final state is \ipf) = a\ k a\_ k (a_ K ) A ' _2 |vac). 
The change of the interaction energy between and 
\tpf) is AEiat = E intJ - E int ,i = <±vU + - - 2vU- , where 
v = N/Nq is the average filling per DW. The energy 
conservation condition reads AE int + AEk in — 0, i.e., 



2e_(K) - e+(ki) - e+(k 2 ) = 4j/U+_ - 2^U_ 



(7) 



By employing the Fermi Golden Rule and the tight- 
binding condition in the z-direction, in which the system 



is tightly confined, we arrive at the expression for the 
transition rate per DW: 



1 / 47ri/ft 2 a i 



J d 2 r 



\I>* 



,K 



kl,k; 



P ( £+ (k 1 ))-l+p( £+ (k 2 ))- 



Here k(r) = e lk r u^ (r) are the Bloch wave-functions 
for the two bands (u^(r) being their corresponding 
lattice-periodic parts), p(e) — J2keB z — £ +(k)) is 
the density-of-states for the '+'-band, and the prime indi- 
cates that the momentum summation is restricted to the 
pairs (ki,k2) that satisfy both the energy-conservation 
condition in Eq. ([7]) and momentum conservation up to 
a reciprocal lattice vector (ki + k 2 = 2K + G). The 
transition rate is calculated by way of numerical evalua- 
tion of the '+'-band density-of-states p(e) and the spatial 
integral / dPr *+, kl *+. k2 * 2 _, K in E q- ©• The latter is 
based on the eigenvectors obtained in the band-structure 
calculation, that is, coefficients in the Fourier expansion 



iG-: 



over reciprocal-lattice vectors. 



To quantify the stability of bosons in the excited band, 
it is pertinent to compare it with the timescale corre- 
sponding to the hopping amplitude th = \Hq\ 
where (j?^ 1-1 are the Wannier functions correspond- 
ing to a pair of adjacent (in the x-direction) DWs. For 
V /E R = -1.0, r = 0.08 we obtain t h = -0.1007 E R . 
With the choice of 87 Rb atoms and A = 810nm (E R w hx 
3.5 kHz), the hopping timescale is h/\th\ ~ 0.45 ms. The 
logarithm of the dimensionless lifetime T — \th\/{fiw), 
depicted with varying filling factor in Fig. [H exhibits 
non-monotonous dependence on the filling. The salient 
feature is that the lifetime is very long for fillings smaller 
than some critical value and also for fillings larger than 
some higher critical value. This is easy to understand 
from the energy-conservation requirement: namely, the 
change in band-energy in the decay process of interest is 
bounded both from above (by 2(e™ m — e™ ax )) and from 
below (by 2(e™ m — e™ m )), while the change in the inter- 
action energy is linear in the filling factor (recall Eq. ([7])). 

The predicted superfluid phase can be experimen- 
tally identified from the time-of-flight density distribu- 
tion (n(r)) t cx E G l*-( k )l 2 52 ( k K G )> where 
k = mbr/(Ht) and $-(k) is the Fourier transform of the 
Wannier function $_(x,y). The resulting density dis- 
tribution, displayed in Fig. [5l has Bragg peaks at the 
condensate wave- vector K and some of the wave- vectors 
related to it by a translation through a reciprocal-lattice 
vector. Importantly, this distribution is asymmetric with 
respect to both k = K and k = 0. This is a combined ef- 
fect of the odd parity of the Wannier function $_(x,y) 
and its extended character in real space (that is, localized 
in momentum space) in the regime of interest. 

The resulting many-body state |\&) spontaneously 
breaks the time-reversal and lattice translational symme- 
tries, with circulating bond currents 17 1 Jij oc (i&lajaj — 
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bosons in the first excited Bloch band of a double-well 
optical lattice. The lifetime of this metastable state is 
shown to be several orders of magnitudes longer than 
the typical tunnelling time. Experimental realization us- 
ing, for example, stimulated Raman transitions 11( or a 
similar method, is clearly called for. 
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FIG. 4: Boson lifetime in the first excited band for different 
filling factors v. Values of parameters are indicated. 
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FIG. 5: Prediction of the density distribution in a time-of- 
flight experiment, for Vq/E r — —1.0, r = 0.08, and conden- 
sate wave-vector K = (0.5, 0.5)fc L . 

ajai|\I / ) oc sin(K • Ry) between sites i and j with rela- 
tive position vector Ry . Unlike the supersolid phase Q , 
which has the superfluid and density-wave orders, our 
state besides the superfluid order supports a current or- 
der. This state does not carry a net current, since the 
group velocity corresponding to a band-minimum equals 
zero; yet, it can potentially have anomalous transport 
properties resulting from breaking the time-reversal sym- 
metry (i.e., violation of the Onsager reciprocity rela- 
tions) [181 ] . 

In summary, we have found that a long-lived non- 
equilibrium superfluid state with broken time-reversal 
and translational symmetries can be realized with cold 
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